THE INTEGRALS IN GRADSHTEYN AND RYZHIK. 
PART 10: THE DIGAMMA FUNCTION 



LUIS A. MEDINA AND VICTOR H. MOLL 



Abstract. The table of Gradshtcyn and Rhyzik contains some integrals that 
can be expressed in terms of the digamma function = dr" logr(as). In 

this note we present some of these evaluations. 



1. Introduction 

The table of integrals [2] contains a large variety of definite integrals that involve 
the digamma function 

d T'(x) 

(1.1) ^ ) = _io g r( a; ) = I ^. 

Here T(x) is the gamma function denned by 

poo 

(1.2) P» = / t x ~ 1 e~ t dt. 

Many of the analytic properties can be derived from those of r(x). The next 
theorem represents a collection of the important properties of r(x) that are used 
in the current paper. The reader will find in [1] detailed proofs. 

Theorem 1.1. The gamma function satisfies: 

a) the functional equation 

(1.3) T(x + 1) = xT(x). 

b) For n G N, the interpolation formula T(n) = (n — 1)!. 

c) The Euler constant 7, defined by 

- 1 

(1.4) 7 = lim \ Inn, 

n — >oo ^ — ' h 

k=l 

is also given by 7 = — T'(l). This appears as the special case a = 1 of formula 
4.331.1: 

7 + In a 



(1.5) 



11 



This was established in [3]. The change of variables t = ax shows that the case 
a = 1 is equivalent to the general case. This is an instance of a fake parameter. 
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d) The infinite product representation 



(1.6) 



r(*>^ ft [(' + !)" 



x/k 



k=l 



is valid for x e C away from the poles at x = 0, — 1, —2, . . 
e) For neNwe have 



(1.7) 

and 

(1.8) 



V 2) 2 2n n\ 



f) For ieC,i^Zwe have the reflection rule 



(1.9) 



T(x)T(l -x) = 



sin 7rx 



Several properties of the digamma function tl>(x) follow directly from the gamma 
function. 

Theorem 1.2. The digamma function 4>(x) satisfies 
a) the functional equation 



(1.10) 

b) For n £ N, we have 
(1.11) 



i>(x + i) = tp{x) + -. 

x 



n-1 1 

V>(n) = + 



fe=i 



In particular, -0(1) = —7. 

c) For ifC away from x = 0, — 1, —2, • • • we have 



(1.12) 



1 00 1 

■00*0 = -7 1- X 7- r- 

a; f-^ fc(a; + k) 

k—1 



fc=0 



x+k k+1 



d) The derivative of ip is given by 

00 1 

(1-13) ^(*) = E77X 



Wx + fc) 2 ' 



k=0 



In particular, ^'(1) = 7r 2 /6. 
e) For neNwe have 



(1.14) 

In particular, 
(1.15) 



n 

^±^ = -7-2^2 + 2^^-^. 

fe=i 

0(i) = -7-21n2. 
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f) For x € C, x £ Z we have the reflection rule 
(1.16) ip(l — x) — tj)(x) + 7TCOt 7TX. 

2. A FIRST INTEGRAL REPRESENTATION 

In this section we establish the integral evaluation 3.429. Severeal direct conse- 
quences of this formulas are also described. 

Proposition 2.1. Assume a > 0. Then 

r°° dr 
(2.1) / [e- x -{l + x)- a ] — =#»)■ 

Jo x 



Proof. We begin with the double integral 

poo pS 

(2.2) J J e~ tz dtdz = I ■ a'-;. 

On the other hand, 

ps poo ps 
I I er tz dzdt= I 



(2.3) / / e- tz dzdt= I - = ln.s. 



i Jo Ji 



We conclude that 

(2.4) / dz = lns. 



00 e~ z - e~ sz 



o * 
This evaluation is equivalent to: 

roo -ax _ -bx L 

(2.5) / dx = ln~, 

Jo x a 

that appears as formula 3.434.2 in [2J. The reader will find a proof in [3J. 
We now establish the result: start with 

s a — 1 l 



r'(o) = / e- s s a ~ l lnsds 



e- s s a ~ x I dzds 

i) Jo 

e~ z j°° s-^-ds- [°° s a - x e-< 1+z Us] — 
10 \ Jo Jo 

This formula can be rewritten as 

dz 



r'( ffl )=r(a) / (e- z ~(i + z)- a ) 

Jo 

This establishes JO- □ 



Example 2.2. The special case a — 1 yields 

1 x 

1 + x J x 



(2.6) / le-*--±-) ^=- 7 . 



This appears as 3.435.3. 

Example 2.3. The change of variables w = — In x gives the value of 4.275.2: 



(2.7) 



l r 



n 



x 



1 



1 - In: 



dx [°° r _o„ s_„n dw 



e- w - (l + w)-" — = -</>(<?). 



In a: Jo w 
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Example 2.4. The change of variables t = l/(x + 1) in (1270 yields 3.471.14: 

r 1 e (l-i/t) _ 
/„ *(1 - 1) 



f 1 e (l-l/*) _ + a 

(2.8) / : - . dt = 



Example 2.5. The result of Example 12.21 can be used to prove 3.435.4: 

< 2 - 9 > r('--iT=)f- h i-''- 

Indeed, the change of variables t — bx yields from (|2.2p the identity 



-6* 1 \ _ r ( -t 1 n dt 



e - — = / e - 



1 + ax ) x Jq \ 1 + at/b J t 



1 + at/fey t 

Formula (|2.5|) shows the first integral is In j and the value of the second one comes 
from (|2~2"|) . 



_ r q\ dx p - q 
e — = 7 



Example 2.6. The evaluation 3.476.2: 

f°° / 

(2.10) ^ (e-* v 

comes directly from f|2 . lj) . Indeed, the change of variables u = x p yields 



Now write 



i=~ r r e - - - j-i *i + 1 r r -J- - ^. 

P Jq \ l + Uj U pj \1 + U ) U 

The first integral is —7 by (|2.6[) and the change of variables v — u q / p gives 



7 1 r 00 / 1 _„\ du 1 z" 00 v - 

11 . 



p gj / s qj v(1 + v)(1 + vp/i) 

Split the last integral from [0, 1] to [l,oo) and use the change of variables x^l/x 
in the second part to check that the whole integral vanishes. Formula (|2.10p has 
been established. 



Example 2.7. Formula 3.463: 

\ 7 

x ~ 2 

corresponds to the choice p = 2 and q = 1 in (|2.10|) . 



(2.11) 
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Example 2.8. Formula 3.469.2: 

AOO 

(2.12) / (e-^ 



Example 2.9. Formula 3.469.3: 
(2.13) 



dx 37 



corresponds to the choice p = 4 and q = 1 in (|2. 10|) . 



dx 7 
x 4 



corresponds to the choice p = 4 and q = 2 in (|2.10[) . 



Example 2.10. Formula 3.475.3: 
(2.14) / (e-* 



= (l-2-") 7 



corresponds to the choice p = 2 n and q = 1 in (|2.10|l . 

The case p = g in (I2.10|) is now modified to include a paramter. 
Proposition 2.11. Let a, b, p e K + . Then 3.476.1 in [2] states that 

, olr , f°° r _„^p dac ln6-lna 

(2.15) 



X 



V 



Proof. The change of variables t = ax p gives 

-ax" -bx"] dx 1 

e — e — = — 
L i x p J a 

Introduce the term 1/(1 + <) to obtain 



bt/a\ & 



1 



J = - 

PJo 
P PJo 



e~*- 



1 \ dt 1 



1 + V i Pjo 
b \ ds 



-bt/a _ 



1 \ dt 



l + t t 



b + as J s 

Adding and subtracting the term 1/(1 + s) produces 

< 2 ' M > if (iT— if)* 

PJo V " + as 1 + s / s 
The final result now comes from evaluating the last integral. 

We now present another integral representation of the digamma function. 



Proposition 2.12. The digamma function is given by 



(2.17) 

This expression appears as 3.427.1 in [2] 



00 / -x p-ax 

■ /„ (-- — 



dx. 
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Proof. The representation (|2.1j) is written as 

/>oo —z 

(2.18) ip(a) = lim / — dz - 



dz 



s z(l + z) a ' 



to avoid the singularity at z = 0. The change of variables z = e* — 1 in the second 
integral gives 



(2.19) 

Now observe that 
(2.20) 



iA(a) = lim / dz — I — 

S ^°J S Z Jln(l+5) 1 



ln(l+«5) e _i 



eft 



< 



d.t 

ln(l+5) * 



0, 



as (5 — > 0. This completes the proof. 

Example 2.13. The special case a = 1 in ([2~T7| gives 3.427.2: 



(2.21) 



1 — e a a; 



e x dx = 7. 



□ 



Example 2.14. The change of variables t = e x in (|2.17p produces 4.281.4: 



(2.22) 



1 



rfi = —ip{a). 



/0 v lni 

Example 2.15. The special case a = 1 in (PT2"2]) yields 4.281.1: 



(2.23) 



1 



1 



/0 K hxt l-t 
Proposition 2.16. Let p, q G R. Then 



= 7. 



(2.24) 



p-i 



9 -l 



/ y lnx 
This appears as 4.281.5 in [2]. 

Proo/. Write 



1 - x 



dx = hip — tp(q). 



x*- 1 \ , r 1 / 1 



(2.25) / (^1 + j ) dx = 

J \lrix l-xj Jo 



r q-l 



In x 1 — x 



dx 



1 xP- 1 - 1 
lnx 



dx. 



The first integral is —ip(q) from (|2.22[) and to evaluate the second one, differentiate 
with respect to p, to produce 



(2.26) 



cP- 1 - 1 



, ■ dx — / x p 1 dx = — . 

dpj lnx J p 



1 



The value at p = 1 shows that the constant of integration vanishes. The formula 
(12.24P has been established. □ 
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3. The difference of values of the digamma function 

In this section we establish an integral representation for the difference of values 
of the digamma function. The expression appears as 3.231.5 in [2J. 

Proposition 3.1. Let p, q € K. Then 

(3.1) / : 

Jo i — X 

Proof. Consider first 

(3.2) 1(e) =[ x'P-^l-xy-Ux- [ afl- l (l- xf- 1 dx, 

Jo Jo 

that avoids the apparent singularity at x = 1. The integral /(e) can be expressed 
in terms of the beta function 

(3.3) B(a,b)= f x a - 1 (l-x) b - 1 dx 

Jo 

as /(e) = B(p, e) — B(q, e), and using the relation 

(3.4) B ^ = Wt§ 

T(a + b) 

we obtain 

r(p) r( g ) 

T{p + e) T(q + e) 



(3.5) /(e)=T(e) 
Now use r(l + e) = eT(e) to write 

(**\ rn rnM ^lM^lk±A_L_ r( g )-r( g + e) l \ 

3.6) J(e) = r(l + e) — — — , 

V e T(P + e) e r (9 + e)/ 

and obtain (|3.1[) by letting e — > 0. □ 



Example 3.2. The special value ip(l) — —7 produces 

(3.7) / — dx = ~f + 1 p(q). 

Jo i — x 

This appears as 3.265 in [2j. 

Example 3.3. A second special value appears in 3.268.2: 

(3.8) f — —x b ~ l dx = i/>(a + b)-il>(b). 
Jo 1 ~x 

It is obtained from (|3.ip by choosing p — b and q = a + b. 
Example 3.4. Now let q = 1 — p in (|3.ip to produce 

f 1 x p ~ l — x~ p 

(3.9) / dx = ipCl — p) — ip(p) = c °t ""p. 

Jn 1 — X 



This appears as 3.231.1 in [2j. 
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Example 3.5. The special case p — a + 1 and q = 1 — a produces 

r 1 x a — x ~ a 1 

(3.10) / dx = ^(1 — a) — + a) — ir cot ixa . 

J 1 - x a 

where we have used (I1.10[) and ()1.16p to simplify the result. This is 3.231.3 in [3J. 
Example 3.6. The change of variables x = t a in (|3.1[) produces 

(3.11) rt^- t ^ dt = ^ q) - m ^ 

y ' Jo l-t a a 

Now let p = 1, a — v and q = - and the replace \x by p and v by q to obtain 3.244.3 
in0: 



/ ^f-* -\ (7 + *(| 



1) 



Example 3.7. The special case p — b/a and q = 1 — b/a in (|3. 1 1|) produces 
(3.13) / ^ 1 ~ ^ " - dx = - m - b/a) - Mb/a)) . 



JO 

The result is now simplified using (|1.16[) to produce 

,„ A , f 1 X^ 1 - X*- 1 " 1 , 7T 1Tb 

(3.14) / dx = - cot — . 

J 1 — x a a a 

This is 3.244.2 in [3J. 

Example 3.8. The special case a = 2 in (|3 . 1 1 1) yields 

rl f2ju-l _ f2i/-l i 

(3.15) y o — r — ^ — di = -(#,)-v(M))- 

The choice /x = 1 + p/2 and ^ = 1 — p/2: 

nl x p - x _p _ 1 
/„ 1-x 2 21 X ~2 

The identities ^(x + 1) = tp(x) -\-l/x and ip(l — x) — ip{x) = tt cot7nr produce 

/* X^ — X ^ 7T f UTT \ 1 

(3.17) / —xdx = -coti^-)--. 

1 ' J 1 - x 2 2 \ 2 J p 

This appears as 3.269.1 in [2]. 

Example 3.9. The choice fj, = 2±1 and z/ = ^±1 in (j3~l"5]) gives 3.269.3: 

1 .x a , 1 / , /& + 1\ , f a + l 

o 1 - x ' 



(3.16) / x dx = - (Ml + p/2) p/2)). 



(3.18) / ±_^=~(^(^)-^ 



4. Integrals over a half-line 

In this section we consider integrals over the half-line [0, oo) that can be evaluated 
in terms of the digamma function. 
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Proposition 4.1. Let p, q £ R. Then 

(4-1) [°° ( - jt^-) - = ^(q) - i/j(p). 

This is 3.219 in [2]. Also 

<4 - 2 » r (jttw - jrhr) t -««>-«•»■ 

Proof. Let £ = x/(l — x) in (13. ip . The second form comes from the first by the 
change of variables x i— > l/o;. □ 

Example 4.2. The special case p — 1 yields 

(4-3) r (t^~ - 7T^-r) ~ = (<?) + 7- 

This appears as 3.233 in [2]. 

Example 4.3. The evaluation of 3.235: 

,A A, [°° {I + XT - 1 dx ... . 

(4.4) I { J x)b - = W)-W-a) 

can be established directly from (|4.3|) . Simply write 

(1 + x) a - 1 dx _ f°° ( 1 1 \ da; Z" 00 / 1 1 \ dx 



{l + xf x J \l + x (l + x) b J x J \l + x (l + x) b - a J x 
to obtain the result. 

Some examples of integrals over [0, oo) can be reduced to a pair of integrals over 
[0,1]. 

Proposition 4.4. The formula 3.231.6 of [2] states that 

f°° x p ~ l — x^ 1 

(4.5) / dx — 7r (cot up — cot irq) . 

Jo 1 - % 

Proof. To evaluate this, make the change of variables t = 1/x in the part over 
[1, oo) to produce 

(4.6) / — dx = / dx - / — dx. 

Jo l-x J l-x J l-x 

Now use the result (13.111 to write 



(4.7) / dx = ^(q)-iP(p)-mi-q)-^(l-p)}. 

Jo i — x 

The relation ^(a:) — — x) = —n cot(7rg) yields the result. □ 

5. An exponential scale 

In this section we present the evaluation of certain definite integrals involving 
the exponential function. These are integrals that can be evaluated in terms of the 
digamma function of the parameters involved. 
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Example 5.1. The simplest one is 3.317.2: 

(5-1) [ ( — r- - — - ) dx = Mq) - <ij}(p) 

that comes from (|4.2|) via the change of variables x i— > e~ x . 

Example 5.2. The special case p = 1 and = —7 produces 3.317.1: 

(5 - 2) £(tt^-(TtW)" i = * (,)+7 

Example 5.3. The evaluation of 3.316: 
r°° (1 _|_ e -x\p _ 1 

(5-3) J (l +e -x)g dX = ^(«) " " 

comes directly from (|5.1I) . 

Proposition 5.4. Let p, q £ R. Then 

(5.4) / — L_dt = ^(g)-^(p), 

Jo 1 — e 

This appears as 3.311.7 in 2\. 

Proof. Make the change of variables x = e~* in (|3.ip . □ 

Example 5.5. The evaluation (|5.4[) can also be written as 

r°° e *(i-p) _ e *(i-«) 

(5.5) / — dt = 0(g) - ^(p), 

Jo e — 1 

Example 5.6. The special case p — 1 and q = 1 — ^ is 



DC' 



(5.6) / ^ = ^,(1-2,) -4,(1), 
Jo e — 1 

and using "0(1) = ~7 and 0(1 — v) = + 7r cot wv, yields the form 
f°° l _ e "t 

(5.7) / — di = tpM + 7 + 7r cot -KV, 

Jo e - 1 

as it appears in 3.311.5. 

Example 5.7. Another special case of (|5.4I) is 3.311.6, that corresponds to p = 1: 



(5.8) / —dt = 1>(q)+y. 

Jo 1 - e 

Example 5.8. The evaluation 3.311.11: 

(5.9, r^^-^-M T -A-*'' " 



follows directly from f|5.4|) by the change of variables t — (r — s)x. 
Example 5.9. The evaluation of 3.311.12: 

(5.10) r^dx= 1 u r inc - in M - v, f lnc ~ lna 

J c x — d x In c — In d \ \ m c — In a/ \ In c — In d 

is proved by simply writing the exponentials in natural base. 
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Example 5.10. The formula 3.311.10 had a sign error in the sixth edition of [2] 
: it appears as 

f°° e~ px - e~ qx , vr / pn 

( 5 - n ) / ; , dx = z—, cot 



It should be 



7T / p7T 



(5.12) / , , . dx = cot. 

Jo l-e-(P+i) x p + q \p + q 

The value IpTS]) yields 

(5.13) r e ~ px ~ e q * dx = _l_ ^ / _ v , 

and the trigonometric answer follows from (|1.16p . This has been corrected in the 
current edition of [2]. 

Example 5.11. The evaluation of 3.312.2: 

p°° (\ _ e ~ ax )(l - e~ bx )e~ px 

(5.14) / i '- dx = ip{p+a)+i/){p+b)-i()(p+a+b)-i()(p) 

Jo 1 - e 

follows directly from (|3.ip . Indeed, the change of variables t — e~ x gives 

f 1 t p - l (l-t a -t b +t a+b ) , 

(5.15) 1 = 1 >* 

and now split them as 

(5.16) /= / dt- / dt 

Jo 1-t Jo 1-* 

and use (|3.1[) to conclude. 

6. A SINGULAR EXAMPLE 

The example discussed in this section is 

fOO -fix J 

(6.1) / \ % = cotfa*), 

J-oo b- e x 

that appears as 3.311.8 in [2 . In the case b > this has to be modified in its 
presentation to avoid the singularity x = — In b. The case b < was discussed in 
[I]. In order to reduce the integral to a previous example, we let t = e~ x to obtain 

, „~» x dx f^t^dt 

(6.2) 



DO 



b - e- x ./ n b-t 



The change of variables t = by yields 

, x f°° e-» x dx ,„ , f°° y^~ l dy 
(6.3) / — =lf- 1 ' 



j~ocb-e x J a 1-y 

Now separate the range of integration into [0, 1] and [1, oo). Then make the change 
of variables y = l/z in the second part. This produces 

poo -fix j pi _ -n 

(6.4) / <—J±=b^\ * Z dz. 

J-oc b-e J Q 1 - z 

This last integral has been evaluated as cot (717*) in 
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7. An integral with a fake parameter 
The example considered in this section is 3.234.1: 

(7 - 1} / (t^x~~^~; dx = ^ cotnq - 

We show that the parameter a is fake, in the sense that it can be easily scaled out 
of the formula. The integral is written as lim 1(e) where 

£<?-! x -1 \ 

dx — I — — dx. 



Jo (l-ax) 1 ^ J (a-x) 1 - 

Make the change of variables t = ax in the first integral and x = at in the second 
one to produce 

1(e) = a- f tq " dt - a""- f" ™ 

and then let e — *■ to produce 

i/^?-i rr-q \ _ ( f a t ^ l dt f 1,a t- q dt 

dx = 



1 — ax a — x J W 1 — t J 1 — t I 

Differentiation with respect to the parameter a, shows that the expression in paren- 
thesis is independent of a. It is now evaluated by using a = 1 to obtain 

dx = a~ q [ / — dt 



1 — ax a — xj \./n 1—t 
The evaluation (|3.1|) now yields 

" dx = a-«(^(l-?)-^(g)) 

= a~% cot 7rg. 



1 — ax a — x 



Formula (17.11) has been established. 



8. The derivative of ip 

In a future publication we will discuss the evaluation of definite integrals in terms 
of the polygamma function 

( d \ n 

(8.1) Poly Gamma [n, x] := ( — J ip(x). 

In this section, we simply describe some integrals in [21 that comes from direct 
differentiation of the examples described above. 



Example 8.1. Differentiating (|3.1|) with respect to the parameter p produces 
4.251.4: 

(8.2) lnX dx = -iP'(p). 
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Example 8.2. The change of variables x — t q in Ij8.2|) . followed by the change of 
parameter p 1— ► ^ yields 4.254.1: 

Example 8.3. Replace g by 2g and p by (j in (|8 . 3|) to produce 

/ o T - w .*=- 3? 0'(i). 

To evaluate this last term, differentiate the logarithm of the identity 

2 2x ~ 1 

(8.5) T(2x) = ^^T(x)T(x + ±), 

V 71 " 

to obtain 

(8.6) 2iP{2x) = 2 In 2 + i){x) + 4>{x + ±). 
One more differentitation produces 

(8.7) 4ip'(2x)=ip'{x)+<i/j'(x+±). 
The value a; = i gives 

(8.8) ^(|)=3^(1) = ^. 
Therefore we obtain 4.254.6: 

(8.9) / — dx = - — . 

V ' J l-x 2 i 8g 2 

Example 8.4. Differentiating (|3.12p n-times with respect to the parameter p pro- 
duces 4.271.15: 

(8.10) f\» x *±te = _ 1 



9. A FAMILY OF LOGARITHMIC INTEGRALS 

Several of the integrals appearing in [2] are particular examples of the family 
evaluated in the next proposition. 

Proposition 9.1. Let a, b £ R+. Then 

(9.1) f x^il-xf- 1 lnxdx= J ^£^ Wa)-t/j(a + b)). 
Jo T{a + b) 

Proof. Differentiate the identity 

(9.2) [\^(l-^d X = ^§ 
Jo T(a + b) 

with respect to the parameter a and recall that T'(x) = ip(x)T(x). □ 



The next corollary appears as 4.253.1 in [2]. 
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Corollary 9.2. Let a, 6, c G K+. Then 

<"> r^-^ lta **-?^ra(*(i)-*(i + '))- 

Proof. Let t = x c in the integral (|9J]l . □ 

Example 9.3. The formula in the previous corollary also appears as 4.256 in the 
form 

(9.4) f in ( -) ^ ldx = \ B (a, ™) L f £±i^ - * (a 

Example 9.4. The integral 



(9.5) 



f ^ = dx= f 1 x 2n (l-x 2 )-^ 2 Inxdx 
A vl-x 2 Jo 



that appears as 4.241.1 in [2J, corresponds to a = 2n + 1, 6 = i and c = 2 ii 
Therefore 

/^^Im , r(n + i)r(|) r(/ 1N1 

(9.6) I -== dx = [*(n + |) - *(n + I)] . 

Using dUTJ), (PH) and ([134]) yields 

V 7 Jo ^/T^ 2 2 »+i ^ fc J 

This is 4.241.1. 

Example 9.5. The integral in 4.241.2 states that 

( 2„)H / 

Writing the integral as 

(9.9) 1=1 x 2n+1 {l- x 2 )- 1 ' 2 \nxdx 

Jo 

we see that is corresponds to the case a = 2n + 2, b = |, c = 2 in (|9.3|) . Therefore 

(9.10) / = r( 4 n r | n 1 | I | y [#i + 1) - ^(n + §)] . 
Using fL~?J) . (|TTTT|) and (fT7l"4|) yields 

/•l„2n+li„™ 9 2n / 2n+l , > fe \ 

<»■"> l^-^y(- + s^) 

This is equivalent to (|9.8p . 

Example 9.6. The integral 4.241.3 in [2] states that 

( , 12) | ,^ = . : ( J <z^ . _!_ _ taa ) . 
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To evaluate the integral, we write it as 



(9.13) /= / x 2n (l -x 2 ) 1/2 \nxdx 



1 



and we see that is corresponds to the case a = 2n + 1, b = §, c = 2 in 
Therefore 

(9.14) 1= r( ^| r 2 ( ) |) + 2)]. 
Using dUTJ), dUT]) and (fTl"!)) yields 

(9.15) / x 2 Vl-3; 2 lnxdx = - n9 _>"/ ln2+- + V i—i- 

Jo 2 2n + 2 (n + l)\ 2n + 2 ^ k 

This is equivalent to ()9.12[) . 

Example 9.7. The integral 4.241.4 in [2] states that 

(9.16) p-^W, . ^ (ln 2 + g -5-L. 
To evaluate the integral, we write it as 

(9.17) /= / z 2n+1 (l-a; 2 ) 1/2 In a: da: 

Jo 

and we see that is corresponds to the case a = 2n + 2, b = |, c = 2 in 
Therefore 

(9.18) I = r( "| [#» + 1) - V(n + §)] . 

41 (n + §J 

Using dUTJ), (fTTTj) and (fLM)) yields 
(9.19) 

, o2n+l / i 2n+l 

(n+l)(n + 2)( 2 „"+ 3 ) ^ 2n + 3 £j fe 

This is equivalent to (|9.1G|) . 

Example 9.8. The integral 4.241.5 in [5] states that 



(9.20) / In a; - x 2 ) 2 "' 1 dx = - ^ , n ffi - Mn + 1) + 7 + In 4] 
Jo 4(2n)!! 

To evaluate the integral, we write it as 

r 1 1 

(9.21) /= / (l-x 2 )"~2 \nxdx 



and we see that is corresponds to the case a = 1, b = c = 2in (|9.3p . Therefore 

(9-22) f= r( ff| B + r if ) [^(D-^n+l)]. 

Using HHl), (fl~TT|) and (fl~14|) yields 

»1 1 (-2n 

(9.23) / (l-a; 2 ) n "2 lna;cfo = 

./n 2^ 



2n\ / n \ 
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This is equivalent to (|9.20|) . This integral also appears as 4.246. 
Example 9.9. The case n = in f9~T)) yields 

(9.24) / 1 ^ = -;in2. 



o 2 

This appears as 4.241.7 in [2]. 

Example 9.10. Formula 4.241.8 states that 

„ x f°° Inxdx _ , „ 

(9-25) / ^__^ = l-ln2. 

J\ x z \/x z — 1 

To evaluate this, let t = 1/x to obtain 

(9.26) I=-[ t(l -t 2 )- 1 ' 2 Intdt. 

Jo 

This corresponds to the case <x = 2, 6=i, c = 2in (|9 . 3|) . Therefore 

(9-27) I = — 

and the value 1 — In 2 comes from (jl.lip and (|1.14p . 

Example 9.11. The case n = in (|9.15p produces 

/•i _ 

/ y/1 -x 2 Inxdx = (21n2 + l). 

Jo 8 

This appears as 4.241.9 in [2]. 

Example 9.12. The case n = in (|9.19[) produces 

/ ir VI - x 2 Inxdx = -(31n2 - 4). 
Jo 9 



(9.28) 



(9.29) 



This appears as 4.241.10 in [2]. 
Example 9.13. Entry 4.241.11 states that 

(9-30) fJ^L= = -^T 

Jo \/x(l-x z ) 

To evaluate the integral, write it as 

(9.31) I=( X - 1/2 (l-x 2 )- 1/2 \nxdx 



o 



and this corresponds to the case a = |, 6 = i, c = 2 in (|9.3[) . Therefore 

r(|)r(§ ) 

4T(f) 

The stated form comes from using (|1.9p and (| 1 . 16|) . 



(9-32) '=^^(1)-^!)] 
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Example 9.14. The identity 



(9.33) / - 1 ==dx = --\n2 



1 x In a; 7r 



o vl - ■'■ 



appears as 4.243 in [2]. To evaluate it, we write it as 
(9.34) I=( x(l - x A y 1/2 Inxdx 



that corresponds to a = 2, b = |, c = 4 in (|9.3p . Therefore, 
(9.35) /= _Lr 2 (i)[^(i)-^l)]. 

The values V>(1) — ~7 ancl ^(5) = — 7 — 2 In 2 gives the result. 
Example 9.15. The verification of 4.244.1: 

is achieved by using (|9.3p with a = | , 6 = ^ and c = 2 to obtain 

(9-37) j = E1^(!)_ 

Using (|1.9[) and (|1 . produces the stated result. 

Example 9.16. The usual application of (|9.3p shows that 4.244.2 is 



(9.40) T(3x) = _ r(g)r(a+|)r(a;+|) 



<M8) /„ Trip 7)1 

where we have used r(|)r(|) = 2tt/V3. It remains to evaluate The identity 

(|1.16[) gives 

(9.39) ^(i)_^(|)=_^. 

To obtain a second relation among these quantities, we start from the identity 

g3x-l/2 

~2tT 

that follows directly from (|1.6[) . and differentiate logarithmically to obtain 

(9.41) il>(3x) = In 3 + - (^(x) + + \) + 4>(x + §) . 
The special case a: = | yields 

(9.42) +^(|) =-27-31n3. 
We conclude that 

(9.43) ^ ( i) = _ 7 _3 ln3 __^ 
and 

(9.44) V , ( |) = _ 7 _3 ln3 + _^. 
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f 1 \nxdx tt ( it \ 

Jo = ~ { Wi) ' 



This gives 
(9.45) 
as stated in 4.244.2. 

Example 9.17. The evaluation of 4.244.3: 

i„i f 1 \nxdx 7T / 7T \ 

proceeds as in the previous example. The integral is identified as 

(9.47) /= ir(|)r(|)[^(|)+ 7 ]. 

The value (|9.44[) gives the rest. 

Example 9.18. The change of variables t = x 4 yields 



(9.48) 



f 1 X p \nxdx = Vn TjE^) r (p+l\ _ (p + 'A 
Jo 16 T(^) Pi 4 J V ' I 4 J 



The last integral is evaluated using (|9.3p with a = 2±1 ; 5 — I and c = 1 to obtain 
/i „d i„ „ j„ /^TYEtl' 

(9.49) 

Jo Vi - ar io i ^ 4 
The special case p = 4n + 1 yields 

(9.50) jT 1 = ^r(n + |) [i>(n + |) - ^(n + 1)] . 

The special case p — 4n + 1 yields 

1 x 4n+1 lnxrfx 0Fl>+±) 

VI -s* = 16»! ^ + " ^ + 1)] • 
Using (fL7| . (fl~TT|) and (fl~14|) yields 4.245.1 in the form 

v y Jo 2 2 "+ 3 ^ £j k J 

The special case p — 4n + 3 yields 

"1 lna . da . 



[V>(ra + 1) - <0(ra + §)] . 



VT^4 16I>.+ §, 
Using (fl~7)l . (fTTTj) and (fTTHE)) yields 4.245.2 in the form 

, n n , ^x^lnxdx 2 2 "- 2 / n 2 ^(-l)A 

Example 9.19. The change of variables t — x 2n produces 

f 1 Inxdx 1 f 1 J_ l 

9.53) / = — / t2«- 1 (i_t)-™ Intdt. 

V ' Jo ^1 - z 2 ™ 4n 2 i v ; 

Then (|9.3p with a = ^- , 6 = 1 — i and c = 1 give the value 
(9 54) /- 1 In a; da; _ r(^) T(l - 1) 1 

1 j 7o vT^^ 4n 2 r(i-i) Ln2nj ni 2 ' JJ 



DIGAMMA FUNCTION 



19 



Using (fl~7) and (fl~T4]l to obtain 

(9 55) l 1 lnxd '' - " L>( ^-^ 7) 

This is 4.247.1 in 2 



o VI - z 2 " 8n 2 sin(^) 



Example 9.20. The change of variables t — x gives 

hixdx 1 Z" 1 J_ i . , A , 

f2n" i (l -t) n hxtdt. 



Using (|9.3p we obtain 

(9.56) T - ^ = r( ^ )r(1 ~^ ) Up (f ) - ^ (1 - f )] 

1 7 io ^»-i(l-a:2) r(l-^j) l ^ K2n ' VK 2n,i 

Proceeding as in the previous example, we obtain 

\nxdx 7r B 



(9.57) 



This is 4.247.2 in [2]. 



In ' 2n/ 



Some integrals in [2] have the form of the Corollary 19.21 after an elementary 
change of variables. 

Example 9.21. Formula 4.293.8 in [2] states that 

r 1 i 

(9.58) / x a ~ 1 ln(l -x)dx = — (ip(a + 1) + 7) . 



This follows directly from (|9.3[) by the change of variables x 1— > 1 — x. The same is 
true for 4.293.13: 



1 

(9.59) / x a - l (l-x) b - 1 \n{l -x)dx =B{a,b)\i>{b)- ip(a + b)]. 
Jo 

Example 9.22. The change of variables t — e~ x gives 

r°° 1 r 1 1 

(9.60) / xe- x {l-e 2x ) n ~2dx = - {l-t 2 ) n -2 \ntdt. 
Jo Jo 

This latter integral is evaluated using (|9.3p as 

(9.61) I = ~ ^^p^ (|) - ^ + 1)) • 



Using (fL7|) and (Tl"JTj) we obtain 

(9.62) jf xe- x (l-e- 2 *) n -zdx = ^ 

This appears as 3.457.1 in [2]. 



v fc=i / 
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10. An announcement 

There are many integrals in [2] that contain the term 1 + x in the denominator, 
instead of the term 1 — x seen, for instance, in Section [3J The evaluation of these 
integrals can be obtained using the incomplete beta function, defined by 

" l t x - x dt 



(10.1) P(x). 

Jo 1 "+" x 

as it appears in 8.371.2. This function is related to the digamma function by the 
identity 



(10.2) p{x) = \ 



These evaluations will be reported in [5]. 

11. One more family 
We conclude this collection with a two-parameter family of integrals. 
Proposition 11.1. Let a, b £ R+. Then 

1 \ dx 7 



(n.i) ^ (, 

independently of b. 
Proof. Write the integral as 

(11.2) / 
Jo 



a 



\ dx 

)t 4 









1 + X b J X 



dx 



The first integral is (a — b)j/ab according to (|2 . 10[) . The change of variables t = x b 
converts the second one into 



<1L3) ik [ e -i +t j t - v 

according to (12. 6|) . The formula has been established. □ 
Example 11.2. The case a = 2" and b = 2 n+1 gives 3.475.1: 

(11.4) / (exp(-a; 2 







1 \ dx 7 



l + x 2n+1 ) x ' 2 r ' 
Example 11.3. The case a = 2" and b = 2 gives 3.475.2: 

(11.5) / (exp(-x 2 



o 



2"\ 1 \ dx 7 



l + x 2 J x 2 
Example 11.4. The case a = 2 and 6 = 2 gives 3.467: 

(11.6) ' ' " ' ^ ' /r 7 



o 



1 + x 2 J x 2 



Example 11.5. Finally, the change of variables t — ax yields 

T(--TTW)f = f^-f(--T^)?- 
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The first integral is In | according to (|2.5|) , the second one is —7. This gives the 
evaluation of 3.442.3: 
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